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DECOUPLING ELASTIC WAVES AND ITS APPLICATIONS 


HONGYU LIU AND JINGNI XIAO 


Abstract. In this paper, we consider time-harmonic elastic wave scat¬ 
tering governed by the Lame system. It is known that the elastic wave 
field can be decomposed into the shear and compressional parts, namely, 
the pressure and shear waves that are generally coexisting, but propa¬ 
gating at different speeds. We consider the third or fourth kind scatterer 
and derive two geometric conditions, respectively, related to the mean 
and Gaussian curvatures of the boundary surface of the scatterer that 
can ensure the decoupling of the shear and pressure waves. Then we 
apply the decoupling results to the uniqueness and stability analysis for 
inverse elastic scattering problems in determining polyhedral scatterers 
by a minimal number of far-field measurements. 
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1. Introduction 

In this paper, we are mainly concerned with the time-harmonic elastic 
wave scattering due to an incident plane wave and an impenetrable scatterer. 
Let D C be a bounded domain with a Lipschitz boundary dD and a 
connected complement M.^\D. D represents an impenetrable elastic scatterer 
embedded in an infinite isotropic and homogenous elastic medium in 
that is characterized by the Lame constants A and n satisfying /r > 0 and 
3A -|- 2/r > 0. Consider a time-harmonic elastic plane wave x = 

I 
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G (with the time variation of the form being factorized out, 
where u G M_|_ denotes the angular wavenumber) impinging on the scatterer 
D. The propagation of the elastic plane wave will be interrupted and this 
leads to the so-called scattering. The elastic scattering is governed by the 
Lame system (or the reduced Navier equation) (cf. [UlllUjl. 

{A*+uj‘^)U = 0 in R\D, (1.1) 


where U{x) = G denotes the total displacement field. In (jl.ip . 

the operator A* is defined by 
3 




1=1 


dxi 




dUi dUi 


dxi dxj 

d 


d 




=^JiAUJ + {\ + lJ)—\V■U] 


( 1 . 2 ) 


d 


— ~ A V A U\j -|- (A -|- 2/i)^—[V • t/], j — 1,2, 3. 


Throughout, the incident/impinging elastic plane wave is of the following 
general form 


17*"(x) = U^'^{x;d,d-^ ,ap,as,uj)apde 


iknd-x 


7_L 

“h (XsCi € 


iksd-x 


:= + UT, d, d^ G Up, as G C, 


(1.3) 


where d G := {x G ||a;|| = 1}, is the incident/impinging direction; 
d-*- G satisfying d • d-*- = 0 is the polarization direction; and ks := uijy/Jl, 
kp := ujjy/X + 2fi respectively denote the shear and compressional wavenum¬ 
bers. In what follows, for a given vector d G S^, d-*- shall always denote a 
unit vector that is perpendicular to d. If / 0 while = 0 for t/*" in 
p.3p . then W"' = Up"' := apde^^p^''^ is the plane pressure wave; whereas if 
Op = 0 and as / 0, then = Ul" := asd^e''^p^'‘^ is the plane shear wave. 
The total elastic wave field U in (ini) is composed of the incident plane wave 
17®” in (II.3p and the scattered wave t/®'®. 


U{x) = W"{x) + U^^ix), X G M.^\D. (1.4) 


It is straightforward to verify that the scattered wave field = U — 17®” 
satisfy the Lame system 

{A* + = 0 in R^D. (1.5) 


The scattered wave 17®®® is radiating and required to satisfy the Kupradze 
radiation condition, 

(V A V A 17"®=)(x) A ^ - iksV A 17"®=(x) = 0{\\x\\-^), 

a, (1.6) 

^ • [V(V • 17*®=)](x) - ikpV • U^^{x) = 0(||x||-2), 

||x|| 

as ||x|| —)■ -|-oo, which holds uniformly in all directions x := x/||x|| for 
X G M^\{0}. 

To complete the description of the elastic scattering problem, we next 
complement the system (ll.ll) - (ll.6l) with a suitable boundary condition on 
the exterior boundary of the impenetrable scatterer D. Depending on the 
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physical property of the scatterer, the following four kinds of boundary 
conditions have been widely considered in the literature. The first kind 
boundary condition is given by 

17 = 0 on dD. (1.7) 

The second kind boundary condition is given by 


TU = 0 on dD, (1.8) 

where T is the stress (or traction) operator on dD defined by 

TU :=2ndyU + \uV -U + ^U), (1.9) 

with 1 / € the outward unit normal vector to dD, and dy the boundary 
differential operator defined as 


dyU := [z. • 

The boundary conditions of the third and fourth kinds are, respectively, 
given by 

u • 17 = 0 and u A TU = 0 on dD, ( 1 - 10 ) 

and 

u A 17 = 0 and v ■ TU = 0 on dD. (l-H) 

The elastic system (in])-(lL6D associated with either of the four kinds 
of boundary conditions in (ll.7l) - (ll.lip is well understood with a unique 
solution U e HI^^{M.^\D). In particular, if dD is C^-smooth, then U £ 
C‘^{M.^\D)^ n C^{M.^\D)^. We refer to [5l[9] for the related results. 

It is widely known that similar to (|1.3p . one has the following decompo¬ 
sition 


use ^ JJSC ^ JJSC^ 


:=--^V(V-17"^), 17f :=-^VA(VA17"'^), (1.12) 

kp kg 

where the vector functions Up^ and U^^ are referred to as the pressure (lon¬ 
gitudinal) and shear (transversal) parts of 17®'^, respectively, satisfying 

{A + kl)U^'^ = 0, VAl7p"^ = 0 in M.\D, (1.13) 

(1.14) 


use 

P 


{A + k'^gp!'^ = 0, V-C/f = 0 in M.\D,. 


Moreover, the Kupradze radiation condition (11.611 is equivalent to the Som- 
merfeld radiation condition. 


= 0(||x||-2), (1.15) 

= 0(||x||-2), (1.16) 


uniformly in all directions x as ||x|| —>■ -|-oo, for j = 1, 2,3. In terms of (II.3p 
and (jl.l2p . one has the corresponding decomposition of the total wave field 
as follows. 


U = Up + Us, Up-.= Wp^ + U;‘^, Us:=ur + u;‘^. (1.17) 

Using the decomposition described above, there holds the following as¬ 
ymptotic expansion as ||a;|| —)■ -|-oo, 

1 




( 1 . 18 ) 









4 


HONGYU LIU AND JINGNI XIAO 


and are known as the P-part (longitudinal part) and S-part (transver¬ 
sal part) of the far-field pattern of Moreover, using the decomposition 
(11.121) . and are the far-held patterns of and respectively. 
In this paper, we dehne the total held pattern to be 

Urix)-.= U^ix) + Urix). (1.19) 

It is easy to verify that 

(x) = [U^ (x) ■ x]x and U^{x) = [x AU^{x)] A x. (1-20) 

Though the elastic wave held 17®'^ can be decomposed into its shear and 
compressional parts, the two types of waves are generally coupled together. 
In fact, by using the boundary conditions (ll.7l) - (|l.lll) . one can readily see 
that the coupling occurs on the boundary of the scatterer D. That means, 
in general, the shear wave and the pressure wave coexist, even in the case 
that the incident wave is a pure plane pressure wave or a pure plane shear 
wave. In this article, we shall derive two geometric conditions in terms of the 
mean and Gaussian curvatures of the boundary surface of the scatterer D, 
that can ensure the decoupling of the pressure and shear waves for the case 
with the third or fourth kind boundary condition. Clearly, if the pressure 
and shear waves are decoupled, then the study of the Lame system m can 
be significantly simplified to the investigation of the two vectorial Helmholtz 
systems (I1.13P and ()1.14p . This undoubtedly would find important applica¬ 
tion in relevant study of the elastic wave scattering. Indeed, another main 
motivation for our this study is the inverse elastic scattering problem which 
concerns the determination of the scatterer D by knowledge of the corre¬ 
sponding far-held pattern U^{x), where r = t,p or s. If one introduces an 
operator T" which sends the unknown/inaccessible scatterer D to the corre¬ 
sponding far-held pattern [/“ associated with the incident/impinging plane 
wave in o, then the inverse scattering problem can be abstractly 
formulated as follows, 

T{D) = U^{x;d,d'^,ap,as,uj,D), x G T = t,poTS, (1-21) 

where we write U^{x;d,d-^,ap,as,uj,D) to specify its dependence on the 
incident plane wave as well as the underlying scatterer. The inverse problem 
(ll.2ip is of signihcant practical importance and has been extensively inves¬ 
tigated in the literature. It is easily seen that the inverse problem (ll.2ip 
is nonlinear and it is known to be ill-posed in the sense of Hadamard; we 
refer to laEllSHZllIZllIB] and the references therein for related results on the 
inverse problem. 

In this paper, we are particularly interested in the case that a single 
far-held measurement is used for the inverse problem (ll.2ip . That is, the 
far-held pattern is collected corresponding to a single incident plane 
wave with d, d-*-, ap, ag and cj being hxed. The inverse problem is formally- 
determined with a single far-held measurement. The inverse elastic scat¬ 
tering problem with a single far-held measurement still remains to be a 
challenging open problem in the literature. To our best knowledge, the only 
existing general uniqueness result in the literature for the inverse problem 
(ll.2ip with a single far-held pattern was derived in a recent work [5] for the 
determination of polyhedral scatterers endowed with the third or fourth kind 
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boundary condition. The argument therein is based a reflection principle 
for the Navier equation associated with the third or fourth kind boundary 
condition, and the path argument developed in for inverse acous¬ 

tic and electromagnetic scattering governed by the Helmholtz and Maxwell 
equations, respectively. In the current study, the polyhedral scatterer is a 
particular case where the pressure and shear waves can be decoupled. Using 
the decoupling result, we are able to derive much more general unique¬ 
ness results which cover the uniqueness results in [5] as two particular cases. 
Moreover, by using the decoupling result, we can advance the relevant study 
much further in deriving optimal stability estimates through the extension 
of the stability results in our recent work m- These are the first stability 
results in the literature in determining elastic obstacles by a single far-held 
measurement. 

The rest of the paper is structured as follows. In Section [2j we derive the 
geometric conditions for the decoupling of elastic waves. Section 3 is devoted 
to a case study when the underlying scatterer is a ball. In Section 4, we 
apply the decoupling result to the study of the uniqueness and stability for 
the inverse elastic scattering problems. The paper is concluded in Section 5 
with some relevant discussion. 


2. Decoupling elastic waves 


In this section, we consider the decoupling of the pressure and shear waves 
for the elastic system (fTTTD-ra associated with, respectively, the third 
kind boundary condition (|l.ini) and the fourth kind boundary condition 
(jl.lip . For notational convenience, we write U (x; [/*"■, [D, III]) = U (x; d, 
ctp, Os, w, [D, III]) and 17(x; [/*”, [D, IV]) = 17(x; d, d-*", a^, a^, a;, [D, IV]) to 
denote the total wave helds corresponding the third and fourth kinds scat- 
terers, respectively. The same notations shall be employed for the scattered 
wave fields and the far-held patterns, as well as the scattering problems asso¬ 
ciated with the hrst and second kinds of boundary conditions. Throughout 
the rest of the paper, we assume that the Lipschitz scatterer D has a piece- 
wise C^’^-smooth boundary. Since the Lame system a is strongly elliptic, 
we know that the solution is locally near a regular piece of 

dD (see [TB]L 

Let r be a regular piece of dD and be represented in the following para¬ 
metric form, 

x(u) = (xl(^^l,^^2),X2(^il,^i2),a;3(^^l,^^2))'^, u = (ui,U 2 )'’" S , (2.1) 


where the superscript T denotes the transpose of a matrix. Denote by 
9 = {9jk)‘j k=i fundamental matrix of differential geometry for T 

with 


dx dx 

9jk * o o 

duj duk 


j,k = 1,2. 


( 2 . 2 ) 


Henceforth, we assume that the parametric form m is such chosen that 


u = u(x) 


1 dx dx 
dui ^ du2 ’ 


\g\ := det( 5 ). 


(2.3) 
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is the outward unit normal vector of dD on F. In what follows, we let 


H{x) = -{ki + K 2 ){x) and K{x) = {KiK 2 )ix), x e F, 


(2.4) 


denote, respectively, the mean and Gaussian curvatures at x G F, where 
Ki and K 2 are the principle curvatures of F at x G F. We also denote by 
Gradr the surface gradient operator on F (cf. [ilfTUj i. and recall that for any 
y? G C^{dD), 


Gradr = 


E 




dip dx 
duj duk ’ 


(2.5) 


where 

■■= l(m)h.ir'■ ( 2 . 6 ) 


2.1. The fourth kind boundary condition. We first treat the scattering 
problem associated with the fourth kind boundary condition and have that 

Theorem 2.1. Let F C dD be a regular piece of the boundary surface of 
the scatterer D. Suppose that there holds 

H{x) = 0 for X G F. (2.7) 

Then for U = a solution to the Lame system (ll.ip - (ll.6p satisfying 

the fourth kind boundary condition on F, 

u AU = 0 and v ■ TU = 0 on F, (2-8) 

if and only if 


V-C/ = 0 


and 


z/ A (V A (V A G)) = 0 on F. (2.9) 


Proof. Let u(x) = (uj(x))j^;^ G be the unit normal vector to x G F 
pointing to the exterior of D. Define 

V[u\{x) ■.= Gradrui(x), GradrU 2 (x), Gradrt' 3 (x) , x G F. (2.10) 


Since 2H{x) = —V • u(x) (cf. [E]), (|2.7I) actually implies that 

3 

= 0 for X G F. 


tr(D[z/](x)) := j^GradrUj(x) 


( 2 . 11 ) 


i=i 


Next, we show that (12.Sp implies (12.9p . By direct calculations, it is easy 
to show that 


Gradr = u;GradrC7j + DjGradru; 

=i'iVUj — V + [/jGradru;, j, / = 1, 2,3. 


( 2 . 12 ) 


Noting that v ■ Gradru; = 0, / = 1, 2,3, then by using the fact that v AU = Q 
on F in ()2.8p . we have 


17 • Gradru; = 0 on F, / = 1,2,3. (2-13) 

and 

Gradr ^ [u A = 0 on F, j = 1,2,3. (2-14) 
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Hence, by virtue of (I2.12I1 . together with straightforward computations, one 
further has that 


0 = 


6 

^ [Gradr(u/[/j-UjC/,) ^ 


1=1 


(2.15) 




(u • - V • t/) + Uj [cradru^j ^ 


i=i 


on r,j = 1,2,3. 


Multiplying both sides of (I2.15p by Vj and then summarizing the resulting 
equalities for j = 1, 2,3, one then arrives at 


C> 

u-d,U = V-U-iu-U)Y[ Gradri 


on r. 


i=i 




(2.16) 


By applying (12.lip to (I2.16p . we clearly have u ■ dyU = V-U on B and as a 
consequence, by using u ■ TU = 0 in (12.8p . there holds 


(2.17) 


0=u-TU = 2/iu -dyU + XV -U 
= (2/i + A)V-t/ on r. 

Finally, it is easy to have that 

A (V A (V A U)) = u A 17 + A fv(V • ?7)'j = 0 on F, (2.18) 

kg kp \ / 

which together with (I2.17p clearly yields ()2.9p . 

We proceed to show that ()2.9p implies ()2.8li . In view of the two conditions 
in (EJl), one readily has 

z. A 17 = A (V A (V A [/)) - A (V(V • 17)) = 0 on F. (2.19) 


1,2 

Ivp 


Furthermore, by using (I2.16p and ()2.1ip . one clearly has v -dyll = V-17 = 0 
on F, which in turn implies that 

u • TU =2/n/ • a^l7 + AV • 17 

= (2/i + A)V-17 = 0 on T. 


The proof is complete. 


□ 


Remark 2.1. It is easily seen from its proof that Theorem 12.11 holds in a 
more general scenario where U is not necessarily a solution to the scattering 
system (fFiD-ifrel). Indeed, Theorem 12.11 still holds, as long as U is an 
solution to the Lame system dni in a neighborhood of F. However, in order 
to have a consistent exposition, we stick to the scattering system ([n])-([L6]) 
for the decoupling results, and the same remark applies to our subsequent 
decoupling study in the present section. 


Remark 2.2. From a geometric point of view, the condition (12.7p means that 
if F is a minimal surface, then the two boundary conditions (12.8p and (12.9p 
are equivalent to each other. 


The following proposition is a direct consequence of Theorem 12.11 
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Proposition 2.1. Let D = [L),IV] he a fourth kind elastic scatterer with a 
piecewise C^’^-smooth boundary such that each of its regular piece satisfies 
the geometric condition (I2.7|) . Let U{[D,IV]) = [/(•; t/*”, [H, IV]) with the 
incident field given by m- 
Set 


Vp = Vp{[D, IV]) := -V • U{[D, IV]), v;^ := -V • IV]). 

Then Vp is a solution to the following problem, 

(A + kp)vp = 0 in 

Vp = i;*” + Vp‘^ in 'Sp\D, 

Vp = 0 on dD, 

77 ^ • — ikpV^^ = 0 (jjxll“^) as Ux]] ^ +oo, 

If!! ^ ^ 

with 

:= -V • 1/7 = -{iapkp)e^^^^-'^. 

In parallel, if we set 

E, = Es{[D,lY]) ■.= V AU{[D,W]), .Ef := V A [/""(p, IV]), 
H, = Hsi[D, IV]) := V A Es/iiks), 
then {Es,Hs) solves the following problem, 

V A Eg — ikgHs = 0 in M.^\D, 

V A Hg + ikgEg = 0 in 

^ Eg = Ef + Ef in M3\;d, 
u A (V A Eg) = 0 on dD, 

^ A (V A Ef) + ikgET = OdlxlJ-^) as Ux]] ^ + 00 , 
||x|| 

with 

Ef := V A t/*” = iagkg {d A 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 

(2.24) 


(2.25) 


(2.26) 


Proof. By using the decompositions in ()1.12p and (I1.17P , it is straightforward 
to show that 

-y -u = -V - Up and V AU = V AUg in R^D. (2.27) 

Hence, by noting (11.131) and (I1.14p . both Vp in (12.211) and Eg in (12.241) satisfy 
the vectorial Helmholtz equation. It is also noted that the Maxwell equations 
in (12.251) for {Eg, Hg) is equivalent to that Eg satishes the vectorial Helmholtz 
equation (cf. H]). The boundary conditions in (I2.2ip and (I2.25p . respectively, 
for Vp and Eg, are immediate consequences of Theorem 12.11 Finally, the 
radiation conditions can be verified by direct calculations. 

The proof is complete. □ 


(|2.22l) is known as the Helmholtz system the describes the time-harmonic 
acoustic scattering problem, where D denotes a sound-soft scatterer and, 
v^jf, Vp^ and Vp are respectively called the incident wave, scattered wave and 
total wave fields. (I2.25P is known as the Maxwell system that describes the 
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time-harmonic electromagnetic scattering problem, where D represents a 
perfect magnetic conducting (PMC) scatterer, and Eg and Hg are 

respectively the incident electric wave, the scattered electric wave, the total 
electric wave and the total magnetic wave helds (see, e.g., a)- In what 
follows, we shall write D = [D, SS] to specify that D is sound-soft scat¬ 
terer corresponding to the Helmholtz system (12.22^ : whereas D = [D, PMC] 
means that D is a PMC scatterer associated with the EM scattering problem 

( 1225 ]). 

It is immediately observed from the definitions (I2.2ip and (11.171) (see also 
dlSZD) that 

UpiW^-,[D,W]) = ^Vvp{[D,W]), 


and similarly. 


Ug{W^; [D, IV]) = ;^V A Eg{[D, IV]). 


(2.28) 

(2.29) 


In what follows, we shall write 


and 


-V • [Z),IV]) =Vp{v;^-, [D,SS]), 

[D,IY]) =^Vvp{v;^-, [D,SS]), 


V A [D, IV]) =Eg{Ef; [D, PMC]), 

Ug{U^^-[D, IV]) =^V A Eg{Ef- [D, PMC]), 


(2.30) 


(2.31) 


to indicate the correspondences between the solution to the elastic scatter¬ 
ing problem and those to the acoustic and the electromagnetic scattering 
problems that are stated in Proposition 12.11 

Due to the symmetry role of the electric and magnetic fields in the 
Maxwell system (|2.25l) . as an alternative to (I2.24p . one can also set 


Eg = Eg{\D,l\]) :=H,([D,IV]), Ef([D,IV]) ■.= Eg-Ef, 
Hg = Hg{[D,l\]) := -Av aE,([D,IV]), 

Ka 


with 


Ef := -—V A (V A t/*”) = -iagkg 

/to 


iksd-x 


(2.32) 


(2.33) 


This would lead to a similar Maxwell system to (|2.25l) for {Eg,Hg), but with 
the homogeneous boundary condition on dD replaced by 


n A Eg = 0 on dD. 


(2.34) 


In this case, [D,IV] is referred to as a perfect electric conducting (PEC) 
obstacle, and we shall write [D, PEC] when it is employed in the EM problem 
(|2.25l) with the boundary condition replaced by (12.341) . Eurthermore, in such 
a case, (|2.31l) should be modihed to be 


Ug(U^^; [D,IV]) = -Es{ET; [D,PEC]). 

/to 


(2.35) 
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Clearly, according to Proposition 12.11 if a fourth kind scatterer D satis¬ 
fies the geometric condition (12.71) . then the corresponding elastic scattering 
system can be decoupled into the two scattering systems (|2.22p and (|2.25p . 
which involves only the pressure and shear waves, respectively, through the 
relations (12.281) and (|2.29l) ; see Theorem 14.11 in Section 0] for a summarizing 
conclusion. This fact would be of significant importance when considering 
the inverse problem of determining the scatterer D from the corresponding 
elastic far-held measurements. Indeed, it can be seen that the information 
of the scatterer is also implicitly contained in the two simpler systems (12.221) 
and (j2.25l) . We shall make extensive use of such a decoupling idea for our 
study in Section [4] for the inverse elastic scattering problems. 

Next, we show that the geometric condition (12.lip on the scatterer [D, IV] 
is also “generically” necessary for the decoupling of the pressure and shear 
waves. Indeed, we have that 


Theorem 2.2. Let D be a fourth kind scatterer such that 

for a.e. x&dD. (2.36) 

Let U be a solution to the Lame system (in])-([L6]). Suppose further that 

in R^\D. (2.37) 

Then (|2.8I) and (12.91) cannot hold true simultaneously on dD. 


Proof. We prove the theorem by an absurdity argument. Assume contrarily 
that both (|2.8I) and (|2.9I) hold on dD for the scattering system (ll.ll) - (ll.6l) . 
First, by using (11.91) . one has by a direct calculation that 

u • TU = 2p,v ■ dyll + AV • U. (2.38) 

Then, by using (12.91) on dD and the fact that V • l/g = 0, one has 

V-U = V ■ {Up + Us) = V ■Up = 0 on dD. (2.39) 

Inserting (I2.39p into (I2.38p . and using (12.81) on dD that u • TU = 0, one 
readily has 

u ■ duU = 0 on dD. (2.40) 

Next, by using ()2.8I) on dD that u A 17 = 0, and a completely similar 
argument to the hrst part of the proof of Theorem 12.11 one can show that 
(I2.16P also holds on dD., namely. 


3 

u-d^U = V -U -{u-U)Y^ [Cradruj 
Clearly, (12.361) implies that 


on dD. 


(2.41) 


tr(P[u](x)) / 0 for a.e. x G dD, (2.42) 

where tr(T>[u]) is defined in (|2.11l) . Combining (|2.42l) . (|2.39l) and (I2.40p and 
(I2.4ip . one can readily verify that 

ly ■ U = 0 on dD, (2.43) 


which together with the fact again that u A t/ = 0 one dD immediately 
implies that 


U = 0 on dD. 
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That is, D is also a first kind scatterer and hence [/(•;[/*"■, [D,IV]) = 
[/(•; [/*"■, [T>,I]) in M.^\D, which is a contradiction to (I2.37p . 

The proof is complete. □ 

Remark 2.3. We would like to remark that physically speaking, (I2.37p means 
that two scatterers possessing the same shape D but of different physical 
properties must produce different scattering effects due an incident plane 
wave 7/*”. Hence, the condition should be a generic one. This is also related 
to the inverse scattering problem of identifying the physical property of an 
unknown/inaccessible scatterer by using a single wave measurement due to 
the incident plane wave 

In Section [3l we shall consider the case that D = Bji, a central ball of 
radius R S M+. In this case, one has that 

H{x) ^ (2.44) 

Hence, the condition (12.361) holds true. By using the method of wave series 
expansion, we show that (12.8|) and (12.91) cannot hold true simultaneously 
on dD. This reinforce our conclusion that the geometric condition (12.111) is 
both sufficient and “generically” necessary for the decoupling of the pressure 
and shear waves of the elastic wave scattering. 

2.2. The third kind boundary condition. In this section, we treat the 
elastic scattering problem associated with the third kind boundary condi¬ 
tion. Our argument shall follow a similar spirit to the fourth kind boundary 
condition case in Section 12.11 but would become a bit more complicated. 
First, we have 

Theorem 2.3. LetT C dD be a regular piece of the boundary of the scatterer 
D. Suppose that there hold 

H{x) = 0 and K{x) = 0 for x G T. (2.45) 

Then for U = a solution to the Lame system (ll.ll) - (|1.6p satisfying 

the third kind boundary condition on T, 


ly ■ U = 0 and 

if and only 

V ^TU = 

on r. 

(2.46) 

u • (V (V • U)) = 0 and 

u A (V A 7/) = 0 

on r. 

(2.47) 


Remark 2.4. The geometric condition (|2.45l) actually implies that T is flat, 
namely v = constant on T. Stating it in the form (I2.45P would give more in¬ 
sights about the geometric condition, especially when treating its “necessity” 
in a certain sense; see ()2.65l) in Theorem 12.41 

Proof. We first show that (I2.46P implies (|2.47l) . By using v - U = 0 in (I2.46p 
on T, one has by straightforward calculations that 

0 =iy A Gradr(u • U) 

3 3 

=1^ A + V /\ 7/jGradrUj 

i=i i=i 

=v A SiyU on r. 


(2.48) 
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where the differential operator is defined by 

:= {v • e C^. (2.49) 

On the other hand, we have 

+ 1^ A{V AU). (2.50) 

Hence, by using u A TU = 0 on F in (|2.46l) and (|2.50l) . one further has 

0 = u A TU = 2^v A — //u A (u A (V A U)) on F, (2-51) 

which together with (I2.48p immediately implies that 

uA(VA 17) = 0 on F. (2.52) 

Furthermore, one has 

u • (V A (V A U)) = Divr(u A (V A U)) =0 on F, (2.53) 

where Divp denotes the surface divergence operator on F (cf. [U Section 3]). 
Thus, one has by (I2.46p and (|2.53l) that 

• (V (V • [/)) = u • [/ - • V A (V A 17) = 0 on F. (2.54) 

kp kg 

Next, we show that (j2.47h implies (I2.46p . First, by using the condition 
u A (V A 17) = 0 on F in ()2.47p . similar to (|2.53p . one readily has 

u-(V A (V A 17)) = 0 on F. (2.55) 

Then using the first relation in (I2.54p . i.e., 

• (V (V • [/)) = u • [/ - • V A (V A 17), (2.56) 

kp kg 

and the fact that u • (V (V • 17)) = 0 on F in ()2.47p , one readily has v-U = 0. 
Hence, there holds that 

u A .^1/1/= u A Gradr(u • 17) = 0 on F. (2-57) 

Finally, by combining (I2.50p . (I2.47P and ()2.57p . one can conclude that 

u A TU = 2/ru A ^^U — //u A (u A (V A U)) =0 on F. (2.58) 

The proof is complete. □ 

In a similar manner to ProDosition l2.ll by applying Theorem 12.31 we have 
the following result. 

Proposition 2.2. Let D = [ZI,HI] he a third kind elastic scatterer with a 
piecewise flat boundary; that is, the geometric condition ()2.45p is fulfilled. 
Let [/([i7,HI]) = [/(•; 1/*”, [77,111]) with the incident field 1/*"' given by (11.31) . 
Set 

Vp = Vp{[D,lll]):=-V-U{[D,m]), v;^:=-V-U^fl[D,lll]). (2.59) 
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Then Vp is a solution to the following problem, 

(A + kp)vp = 0 in Mp\D, 

Vp = + Vp^ in M.^\D, 

V ■ Vvp = 0 on dD, 

77 ^ • — ikpV^^ = C>(||x||“^) as ||x|| ^ +oo, 

IfII ^ ^ 

with v^jf given by (I2.23p . 

In parallel, if we set 

Es = Es{[D, III]) := V A U{[D, III]), Ef := V A W^{[D, III]), 
Hs = Hs{[D,lll]) :=V A Es/{iks), 
then {Es,Hs) solves the following problem, 

A Eg — ikgHs = 0 in 

V A Hg + ikgEg = 0 in M.^\D, 

^ Eg = Ef + Ef in R^\D, 
u A Eg = 0 on dD, 

^ A (V A = 0(11x11-2) ii^ii ^ 

||x|| 

with Ef^ given by (|2.26l) . 


(2.60) 


(2.61) 


(2.62) 


The Helmholtz system (I2.60p is known as the scattering for acoustic waves 
from a sound-hard scatterer. Hence, [D, HI] can be regarded as a sound-hard 
scatterer, which shall be denoted as [L),SH] in this situation. Moreover, in 
parallel to (I2.25p . the Maxwell system (12.621) describes the scattering for 
electromagnetic waves by a PEC scatterer [D,PEC] = [D,HI]. 

In what follows, similar to (I2.30p and (|2.3ip . the correspondences between 
[/(•;[/*"', [D,HI]) and the solutions to the problems (I2.60p and (I2.62p shall 
be written as 


and 


-V ■ -,[0,111]) 

Up{U^^-,[D,m]) 


=Vp{E;^-,[D,SR]), 

Kp 


(2.63) 


V A [D,m]) =Eg{Ef-, [T>,PEC]), 

Ug{U^^-, [D, HI]) =^V A Eg{ET-, [D, PEC]). 


(2.64) 


Next, we show that the geometric condition (12.451) on the scatterer [D, HI] 
is also “necessary” in a certain sense for the decoupling of the pressure and 
shear waves. In fact, we can show 


Theorem 2.4. Let D be a third kind scatterer. Suppose that there holds 

K{x) ^ 0 for a.e. x G dD. (2.65) 
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Let U he a solution to the Lame system (ini-drsD. Suppose further that 

inM.^\D. (2.66) 

Then (j2.46D and (j2.47p eannot hold true simultaneously on dD. 

Proof. We prove the theorem by an absurdity argument. Assume contrarily 
that both ()2.46p and (I2.47p hold on dD for the scattering system (jl.ip - (jl.6l) . 
Then we have by from u A TU = 0 and u A (V A t/) = u A (V A Ug) = 0 on 
dD that 


2/zu A d^U = u A TU — A (u A (V A U)) =0 on dD, (2.67) 
which together with (I2.50p readily implies that 


V/\ = 0 on r, (2.68) 

where is defined in ()2.49l) . Next, by using u • t/ = 0 on dD in (I2.46p . we 
have 


0 =v A Gradr(u • U) 

3 3 

=u A + u A t/jGradrUj 

j=i i=i 

3 

=v A + +u A ^ t/jGradrUj, 
1=1 

which together with (I2.68P further yields that 

3 

V A C/jGradrUj = 0 on 5D. 
1=1 


(2.69) 


(2.70) 


Let {gki)k i=i be the hrst fundamental matrix of the differential geometry 
for r, see (12.2p . Using (12.5p . one can verify directly that 




k,l=l 


gii 


dx 

duo 


duk 


E 

k=l 

duj dx 

dui du2 du2 dui ’ 


ki dvi 9 x 
duk dui 

dx 
dui 

dx 


gi2 


dx 

dk,i -r, - Ok,2 — 

du2 dui 

duj dx 


(2.71) 


i = 1,2,3. 
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Thus, we have 


A C/jGradr 
i=i 

3 


1 


- U.- 


'duidu2 


du2 dui 

1 


du , 

^ ' 
‘ k=i 


k dx 


from which we can directly compute that 

2 

CM 


^ du 

dui du2 


^ g A 

k,i=i '' 


(9u/ 




(2.72) 


1 f ^ dv '\ dx 1 / dv \ dx 

V du2 J dui ^ V ) du2' 

In view of (I2.70p and ()2.72l) . one readily has 

U-^ = 0 on dD, 1 = 1,2. 
dui 

Next, we show that the geometric condition (j2.65p implies that 
dv dv 

—— A -— / 0 for a.e. x € dD. 
dui du2 

Indeed, one can represent the tangential vectors dv/dui, 1 = 1,2 as 

2 


(2.73) 


(2.74) 


(2.75) 


(2.76) 


Clearly, (I2.76P implies that (I2.74p holds true. That means the three vectors 
V, dv/dui and dv/dv 2 are linearly independent a.e. on dD. Hence, by 
virtue of (j2.73l) and v ■ U = 0 on dD in (j2.46p . one immediately has 


U = 0 on dD. 


(2.77) 


That is, D is also a first kind scatterer and hence [/(•;[/*"■, [T*, III]) = 
[/(•; [/*”, [T),I]) in M.^\D, which is a contradiction to (I2.66p . 

The proof is complete. □ 


3. Elastic scattering from a ball 

In this section, we consider the elastic scattering from a ball; that is, we 
assume that D = Br. In this case, one has 




(3.1) 


Hence, both the geometric conditions (12.71) and (12.451) are not satisfied. We 
shall show that the pressure and shear waves cannot be decoupled, which 
reinforce our theoretical results in Section [2l 
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3.1. Spherical wave functions. Our arguments rely essentially on the 
expansion method with the spherical wave functions. We first collect some 
preliminary knowledge for the subsequent use, and refer to mES] for more 
related results. 

Let jn(t) and hn{t), n G N U {0} and t G M, be the spherical Bessel 
function and the spherical Hankel function of the first kind, respectively. Let 
{l^(x); m = —n ,..., n} be the orthonormal system of spherical harmonics 
of order n G N. The tangential fields on the unit sphere 

— iGrady^(x) and — ix A Grad 

Vn(n + 1) " ^ ^ ^nin +1) ” ^ 

for m = —n,..., n, n = 1, 2,..., are called vectorial spherical harmonics and 
form a complete orthonormal system in L|(S^) (cf. [1]). 

For n G N U {0}, m = —n,..., n, and / = j or / = h, we introduce the 
scalar functions 

u^[f-,k]{x):=fn{k\\x\\)Y:^{x), (3.2) 

and the vector-valued functions 

iVC[/;A:](x) := V A (x/„(fc||x||)y™(x)) . (3.3) 

It is known that (cf. [3]), u^[j;kp] (resp. ul^[h;kp]) is an entire (resp., a 
radiating) solution of the Helmholtz equation 

(A + kl)v = 0, (3.4) 

and M™[j; kg] and V AM™[j] (resp. M^[h; kg] and V AM^[h; kg]) are entire 
(resp., radiating) solutions of the Maxwell equation 

{A + k‘^^)E = 0, V-E = 0. (3.5) 

By direct calculations, one has for n G N, m = —n,..., n and f = j or h 
that 

VM™[/;fc](x) = kf^{k\\x\\)Y^{x)x + -^fn{k\\x\\)Gra.dY^{x), (3.6) 


M-[/; A:](x) = /„(fc||x||)GradA x, (3.7) 

V A[M::^[f;k]{x)) =nin + l)^^:^^^^Y:\x)x + ^^^ 

(3.8) 


where the function and fn is defined by 

fn{t) ■■= fnit) +tf^{t). 


In what follows, when there is no ambiguity, we omit the letter k from the 
notations [/; k] and M™[/; k] and simply write them as u^[f] and M^[/], 
respectively. 

Recall from [3] that any radiating solutions to the Helmholtz equation 
(13.4p and Ef^ to the Maxwell equations (13.5D for ||x|| > R > 0 have the 
following expansions 

oo n 

n=0m=—n 


(3.9) 
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and 

oo n 

= E E (3-10) 

72=1 rn=—n 

which (along with their respective derivatives) converge uniformly on com¬ 
pact subsets of 11 a: 11 > R. Correspondingly, the far field patterns for the 
radiating solutions and of the expansions (j3.9|) and (l3.1Up . respec¬ 
tively, are given by 


OO 


kn ^ 

R 72=0 m=— 


(3.11) 


and 

E^ix) 


- OO 72 - 

^E E ^(^^s5(rGradFr(^)-C^AGrady™(x)). (3.12) 

^ 72=1 m=—n 


Let 


^>(x,y) := 


I gifclk-yll 


dvr llx — y\\ 

be the fundamental solution to the PDO — (A -|- A:^). For any vector p S 
we have the following expansion (see for instance, [H Theorem 6.29]): 


p^{z,y) =ikJ2 E 


1 


—J ^ n(n-l-l) 

= 1 rn=—n ^ ' 

. oo n - 

+ iEE 


p ■ M^[j\{y) M’^[h]{z) 


k ^ ^ n(n + 1) 

72=1 m=—n ^ ' 


p.VAM™[j](!,)VAAC['>l{2) (3.13) 


+ i:E E P-V<>[j](!,)VK™[A](=), 


n=0 m=—n 


which is uniformly convergent on any compact subset of ||z|| > 
We derive the following lemma for the later use. 


Lemma 3.1. Let d G and p G be two vectors such thatp-d = 0. Then 
one has 


pe 


ikx-d 


oo n n _ 

E E n(n + lf ' ^ 

72=1 m=—n ^ 

1 ^ ^ y72+l _ 

E n(n + lf Gvad Y^id) V A M^[j](x), 

72=1 m=—n ^ 


(3.14) 


and 

^ OO 72 

= -4vr-j; e^^Y^Vu^[j]ix). (3.15) 

n=0 m=—n 


Proof. The proof basically follows the strategies in [H Chapter 2]. For any 
r G M+, we set z = rx and y = — ||x||d. Notice that 

y™(-d) = (-l)^y™(d) and Grady™(-d) = (-l)"+^Grady™(d). 
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From (I3.6I) - (I3.8I) . we have 


P ■ ^in[3\{y) = {-^Tjn{k\\z\\)p ■ (Grady„"*(h) A h), 

p-VAM™[j](y) = {-lY+^^jn{k\\z\\)p-GT&(\Y^{d), (3.16) 

P-"^Un[j\{y) = (-l)”’^^i|^in(fe||z||)p-Grady™((i). 

Passing to the limit r —)• +oo in the expansion (I3.13h for <l>(z,y), inserting 
(j3.16D , and making use of the relations (I3.1ip and (|3.12p , one can derive that 




Akx-d 


CO n n 

V ^- 

^ ^ n(n + 1) 

n=lm=—n ^ ^ 


p ■ (Grad A d)jn{k\\z\\)GTa.d Y^{z) A z 


An-\-l 


^ oo n 

-T T 

k ^ ^ n(n + 1) 

n=lm=—n ^ ' 


p • Grad Grad (i) 


. oo n 

^ Grady™(d) 


jn{k\\z\\) 




n=lm=—n 


which immediately yields (|3.14p . 

Replacing the vector p by d and applying similar arguments as above, one 
can obtain (13.151) . The series (|3.15p can also be obtained straightforwardly 
by taking gradient on both sides of the expansion (see for instance, [H Ghap- 
ter 2] 

oo n 

n=0m=—n 

The proof is complete. □ 


3.2. Elastic scattering from a third or fourth kind ball. Let C/*” = 
LA” + L/]” be the incident plane elastic wave given by the general form (|1.3I) . 
In the rest of this section, we shall always let u(x) = x be the outward unit 
normal vector to the central ball R/j. 

For n G N U {0} and m = —n,..., n, we define 

Al-(l) := 47rz”+ii™, A^{2) := ■ Gradl^, 


and 

A™(3) := 47ri”d-*- • ^Grady^(d) A . 

It is directly computed from (I3.15P and (13.6p for 17^ = apde^^^'^'^ that 

OO n 

x-UY= -apE E A^Y)UkvM\)yY{x), 

n=0m=—n 

xAUY= -apf^ X; ^;r(l)4^|^^AGradyr(^), 

n=0m=—n P'' 

and that 

oo n 

V • L/- = OpV • (de*"-''-) = apkp kp], (3.18) 

n=0 m=—n 


(3.17) 
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We also have from (13.141) . (|3.7I) and (|3.8I1 for Uf^ = that 

x-UT= -c.f; XI 


72=1 m=—n 
oo 


fcJlxl 


xAUi^= asY^ ( (3)jn(A:s||a;||)Grady™(x) (3.19) 

nin + 1) 

n=l m=—n 


—a. 


E 


—; nin + 1) ^ 

=1 ' ' m=— 


Y 4 l-( 2 )^ 4 M#^AGradyr(^)- 


k„\\x\ 


It is obtained by (13.9p and (13.10^ that, for any radiating solution of 
the Navier equation dEU) outside a certain ball, = Up^ + has the 
following series expansions 


OO n 




n=0 m=—n 
oo n 


(3.20) 


E E (CVAMr[/i;fe,] + CMr[h;fc,]). 

72=1 rn=—n 

Using the series representations in (j3.20p . it is directly computed that 


and 


OO n 


E «V‘-,kp\, 


n=0m=—n 


(3.21) 


and 


x-U^^= kpY Y O'nh'nikpWxlDY^ix), 

72=0 m=—n 
- oo 72 

^ n(n + l)6-/i4fc,||x||)yrC 

'' II 72=1 m= —72 


X 


XI E 


72=0 m=—n 
oo 72 


A Grad y™(x). 


X AUg^ = Y^ E '^™^f2(^s||a^||)Grad y)("(a 


72=1 m=—n 
oo 72 


+ E E C ^''^^PE AGrady„"^(x). 


(3.22) 


(3.23) 


72 = 1 772= —72 " " 

The following lemma shall also be needed. 

Lemma 3.2. Let be defined in (I2.49P with p(x) = x. Then there holds 

SiylJix) = V(u • U){x) — 7 E-U(x) + -;^v{x)iv ■ U){x), X e dBji. 

lull ll^ll 

Proof. It is clear that 

ly. djU = dj{v-U)-U ■ djv, j = 1,2,3. 
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The proof can be concluded by noting for j = 1,2,3 that 


dju{x) = -n—“ TT^ 
X xr 


u. 


□ 


3.2.1. The fourth kind boundary condition. Now we consider the case that 
the ball Bji is a fourth kind scatterer. By using the series expansions in¬ 
troduced earlier, we determine the scattered field = Up^ + in the 
form (13.201) . such that the fourth kind boundary condition (11.111) holds for 
U = on dBji, namely, 

1 /A =-u A and u ■ =-v ■ TW'^ on BBr. 

For the condition u f\U = 0 on BBr., one first sees from (13.171) . (I3.19|) and 
that 

m _ _ ^^(3) jnjksR) /„ „ .X 

^n{n + l)hn{ksR)' 

On the other hand, from Lemma 13.21 and the equation (j2.50l) . it is verified 
that 

u -TU = 2fiu ■ + XV -U 

= 2nn-V{n-U) + XV -U. 

Inserting (13.181) and (|3.17l) into (I3.25p . we then arrive at 
u • TU^p^ = 2/x u • V(u • up + AV • 

oo n 

= ^pY1 (l)fcp(^Jn(A;p||x||) - 2/rj"(A:p||x||))y„”"(x). 

n=0 m=—n 

In the same way, applying (I3.19P to (I3.25p . we have 


OO n 


u • Tl/r = 2/iu.V(u •[/:-) = a. 

n=lm=—n 

where the functions jn, n & N, are given by jn{t) '■= jn{t) — tjp)- Similarly, 
we have from p3.2ip . (I3.22p and (I3.25p that 

u-TT/f = 2/rx-V(x-t/f) = -2/rf; Cn(n + l)^^M^y-(x), 


n=l m=—n 


and 


jSC 


u ■ TWp^ = 2/rx • V (x • Up + AV • U^p 

OO n 

Y C^p(^^n(A:p||x||)-2/i/r"(A:p||x||))y™(x), 


n=0 m=—n 


with the functions hn, n € N, defined by hn{t) := hn{t) — th'p). 

If both p2.8n and (12.91) hold for the scatterer [Bij,IV], then one should 
have n A Up = n A Us = 0 on OBr. As a consequence, it is deduced from 
P3.18D and (I3.23P that a™ = off and bff = 6™ with 


S- := apA^il) 


1 jn{kphT) 
^p ^n(.^pBf 


and 6™ = a. 


A;^(2) 1 UksR) 

n{n + 1) kg hn{ksR) ’ 
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for n G N+ and m = —n... ,n. 

Define Up := and Us := U^ + U^^ with 

oo n 

E E <^<[h-,kp], 

n=l m=—n 
oo n 

n=l m=—n 


(3.26) 


where the constants c!^, m = —n,..., n, n G N, are given by (|3.24|) . Then 
U := Up + Us is a solution to the exterior problem for the Navier equation 
(tni, which satisfies (|2.9p on the boundary dB^. 

By using the property 

jn{t)Kit) - j'n{t)hn{t) = l/t^ 


(cf. mm ), one can directly compute for x G BBr that 


V ■ TUp{x) 


E 


n=l m=—n 


kj KikpR) " ^ 


and that 


* OO 7X i fry, / C\\ 


^ k1hn{ksR)' 

n=lm=—n * ' ' 


Therefore, noting that u • TU = 0 on dB^, in order for the decoupling (I2.9D 
to be true for the elastic scattering problem with [Bji, IV], there should hold 
that 

d,-^ ■ GvadYj^{d) c^p/kp hn{ksR) 


Y^{d) 


ds/k'i hn{kpRy 


m = —n,..., n. 


(3.27) 


for any n G N. Clearly, this cannot be true for all n G N. Hence, (|2.8|) and 
(|2.9p cannot hold simultaneously for the scatterer [H/j,IV]. Moreover, by 
summarizing our earlier discussion, we have the following result. 


Proposition 3.1. For the scattered elastic field of the scattering problem 
(HID-dlSl) corresponding to [H/{,IV], we have 

OO n 

U^y-,[BR,IV]) = U^'^ + J2 E ay{^uy[h-kp]-VAMff[h;ks]), 

n=l m=—n 


where U^^ := U^^ + Uf'^ is given by (I3.26p . and oFff, m = —n,..., n, n G N. 
are constants satisfying 


a. 


hn{tp) - 2n{n + 1) ^”^^^^^^ - 2-^{hn{tp) - h'f{tp)) 


h2 

til e 


4i (as Af^{2) ap A^{l) \ 

klR^ \kl hn{ts) kj hn{tp)) 


m = —n,..., n, n G N, 


with tp = kpR and ts = kgR. 

Particularly, I7(-; [i?/j,IV]) satisfies the decoupling property (|2.9p if and 
only if Oiff = 0 for all m = —n,..., n and n G N. This cannot he true unless 

as = Op = 0. 
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3.2.2. The third kind boundary condition. In this subsection, we consider 
the case that is a third kind scatterer. Similar to our study in Sec¬ 
tion 13.2.11 we determine the scattered held by using the 

series representation (I3.2U|) , such that the total wave U = 17*"' -|- satishes 
(j2.8l) on dBn- Then we show that the decoupling property (12.471) does not 
hold for the scatterer [i?i^,III]. 

In an analogous manner to our argument in Section 13.2.11 we hrst deter¬ 
mine functions tjp := 17®'^ -|- 17*"' and IJs := 17®** — I/*" with 


p 

oo n 

■= Y Cgrad'u™[/i;A;p], 

n=0 m=—n 
oo n 

n=l m=—n 


(3.28) 


such that U = lJp + JJs satishes the decoupling condition (12.471) on dBn- 
It is directly computed that 


OO n 


A{VAUr)=klJ2 Y KK{ksM)GTadY;rix) 


n=l m=—n 
oo n 


+ E E 

n=l m=—n 

and obtained from (|3.14p that 


7c 


■ X A Gradl)("(x). 


OO n 


. A (V A UT) =as y y ^n(3) 7n(fc.|N|) . ^ 

n=lm=—n ^ ^ " " 


oo n 


- a.c 


4m/o\ 

E E , i, ksjn{ks\\x\\)GmdY;r{x). 


(3.29) 


^ ^ n(n -I- 1) 

Thus, by the condition u A (V A Ug) = 0 on dBji, one has 


6 ”* = 


g, 7i;(*(2) jnjksR) 
kg n{n + 1) hn{kgR) 


and = -agYyY_jyMl, 
n{n + l)hn{kgR) 


Moreover, it is observed from (I3.17p and (I3.28p along with the relation v 
ITp = 0 in (|2.47p on SBr that 




jnjkpR) 
kp k!^{kpR') 


So far, we have obtained the total wave held 17 = 17p -|- in (I3.28P which 
satishes the decoupling relation (I2.47P on OBr. Next, we assert that 

U ^U{-,[Br,111]). (3.30) 

Noting that u • 17 = 0 on OBr, we can hrst see from Lemma 13.21 that 

V A%U = u AViv- U) - AU = AU. 

X X 
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Thus, using (|2.50p . one observes that 

v A TU = —2/i -j—^ uAC/ — ^z/A(uA(VA U)) on dBji. 


Define a new operator T on dBn by 

TU := —ij"^ — fii' A (V A U). 

It is easily seen that 

u A TU = u A TU = 0 on OBr. 

Next we establish the representation of v ATU. Notice that u A (V A D) = 0 
on BBr. Thus we have 

V A TU = —-u A U. 

IfII 

Using (13.171) . (13.191) and (13.281) . one can derive on OBr with straightforward 
computations that 


^ OO 

A fUs = - ^ 


R ^ n{n + 1) ^ f^hnits) 

n=l ^ ’ m=-n ^ 


E 7^ElT^AGradyr(^) 


^ OO 


R^n{n + 1) ^ tshnits) 

n=l ^ m=—n * 


E 7lH7TGradyr(^), 


(3.31) 


and that 


ATUp = 2 ifiap-Y^ X AGrady„”"(x). 


A^{1) 


n=0 m=—n R hk y/ 


(3.32) 


It is obvious that v A TU is not identically zero on OBr, which verifies our 
assertion (j3.30l) . 

By summarizing our earlier discussion, we have that 


Proposition 3.2. There hold 

OO n 

u;^i-,[BR,m]) = u;<^ + Y E Pn^<[h-,kp], 

n=lm=—n 


and 

t/f(s[BR,III]) 


^" + E E {7rVAM-[h;fe,] + CM;r[/i;y}, 

n=l m=—n 


where Up^ and U^^ are given in (|3.28p sueh that Up := U^ + and 
Us := t/]"' + Ug'^ satisfying the deeoupling condition (12.471) on OBr, and 
the constants (dff, 7 ™ and C™ are implicitly given by 


tph'n{tp)/3^ - n{n + 1 )/i„(G) 7 ™ = 0, 

2hn{tp)f3'ff + {tshnits) + 2hn{ts)) An 


^ hap A^{1) 
R \kp tph'^itp) 


k‘^ hn{ts) ) ’ 

(3.33) 
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and 

.m ^ ,, A^{3) ^ as 

n{n + 1) t^hn{ts)hn{ts)' 

with tp = kpR and tg = kgR. 

Furthermore, the total field [/(•; [i?/j,III]) satisfies the decoupling property 
(j2.47D on BBr if and only if ffiff = 7 ™ = 7 ™ = 0 for all m = —n,... ,n, and 
n E N. This cannot be true unless Og = Qp = 0. 

Proof. First, one has on dBji that 

00 n 

u-Up = iytjp + n-Y, ^ /3™V<[h;fep] 

n=l m=—n 

00 n 

= E E PnkphnikpR)Y:r{x), 

n=l m=—n 

and 

00 n 

v.Ug = vUg + u-Y, E 7rVAM-[h;A:,] 

n=l m=—n 

=l:»(n+i) x: 

n=l m=—n ' ' 

Hence, the first equality in (I.S.33p can be obtained from v-U = 0 in the third 
kind boundary condition ( 11 . 101 ) on dB^. 

Next, it is observed that 

- 00 n 

^^f{Up-Up)=-2^^-u^Y, E fin^<[h-.kp] 

n=\ m=—n 

= -/iuAE E X A Grad (x), 

n=l m=—n 

and 


uAT{Ug- Ug) 

/ 2 N 00 n 

-/ruAf^ + uAVAjE E {in^ ^ Mff[h;kg] + CXih'^kg]} 

\ II II / m — 1 m — — r) 


c» n 

= -/^E E In (k'ihnikgR)+ 2 

n=lm=—n 
00 n 


n=lm=—n 

h{kgR) 

R^ 


X A Grad Yffi {x) 


E C^#^Grady”(i). 

n=l m=—n 


R 


By using the above facts, the rest of the proof can be concluded by using 
the boundary condition v A TU = 0 in (|1.10l) on BBr, along with the use of 


(non and (lO^ . 


□ 
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4. Inverse elastic scattering problems 


In this part of the paper, we consider the inverse elastic scattering prob¬ 
lems given in ()1.21l) with the elastic scatterer being [L);III] or [D;IV]. We 
call [D;III] and [D;IV] admissible if they satisfy the geometric conditions 
dZIl) and (|2.45p . respectively. The crucial ingredient is the following decou¬ 
pling theorem which can be readily summarized from our study in Section [2j 


Theorem 4.1. Let D be an admissible scatterer and let U be the solution 
to the Lame system (jl.ll) - (|1.6l) . Define 


Vp := —V • U and £’« := V A U. 


(4.1) 


Then there hold 

U = Us + Up, 

Moreover, set 


Us = A Es 


and 


~ U2^'^P- 

Up 


(4.2) 


:= -V • and Ef := V A U^^. (4.3) 

Then Vp satisfies an acoustic scattering system (cf. (I2.22p or ()2.60p ), and 
Es satisfies an electromagnetic scattering system (cf. (I2.25P or (I2.62p ). The 
correspondences are given as follows, 


vp{v^- [Z1,C(K)]) = -V • U{W^; [I1,K]), 
Es{ET-, [D,B(K)]) = V A UiW^-, [D,K]), 
where the notations C(K) and B(K) are defined by 

C(K) = SH and B(K) = PEC, if K = III; 
C(K) = SS and B(K) = PMC, if K = IV. 


(4.4) 


(4.5) 


The next proposition establishes the correspondences between the far- 
fields patterns of the decomposed wave fields in Theorem 14.11 which can be 
obtained by direct calculations. 


Proposition 4.1. Let U, Us, Up and Vp, Es be the scattering wave fields given 
in Theorem \f. 1\ Let Uf^,Uf° and U^ be given in (jl.lQp . namely, 

ur{x)-.= u^{x) + ur{x). (4.6) 

Let and E^ be the far-field patterns, respectively, of Vp and Es that are 
read off from the following asymptotic expansions as ||x|| ^ -|-oo, 

Vp{x) = v^jfiix) + v^{x) + O > (4-7) 

and 

l/k 11 IT 11 / 1 \ 

Es{x) = ET{x) + ^-^ET{x) + O ^ . (4.8) 

11*11 Vll*ll / 

Then there hold 

U^ix) = f *“(*)*, U^ix) = Ua Erix), (4.9) 

Kp Kg 

and 

{x) =—ikpX ■ U^(x), Ef^{x) = iksX f\Uf°{x). (4-10) 




















26 


HONGYU LIU AND JINGNI XIAO 


Clearly, by using the correspondences given in Theorem 14.11 and Propo¬ 
sition the relevant results for the inverse acoustic and electromagnetic 
scattering problems can be extended to the inverse elastic scattering prob¬ 
lems with admissible scatterers of the third and fourth kinds. We are par¬ 
ticularly interested in the inverse problem of determining D in (11.2111 by 
knowledge of a single far-held measurement; that is, U^{x\ C/*"'), t = t,p or 
s, is collected corresponding to a hxed incident plane wave 17*"'. There is a 
widespread belief that the uniqueness can be established for such an inverse 
problem. However, it still remains to be challengingly open. As discussed 
in Section [H in a recent work [5] , the authors established two uniqueness 
results in determining polyhedral scatterer of the third and fourth kinds. 
By using the decoupling results in Theorem 14.11 and Proposition 14.11 it can 
be shown that much more general uniqueness results hold for the inverse 
elastic scattering problems with polyhedral scatterers. Moreover, we can 
step much further in driving optimal stability estimates for the inverse elas¬ 
tic problem through the extension of the recent stability results in m for 
the inverse acoustic scattering problems. It is emphasized that depending 
on the presence of screen-components or not, we shall make use of a single 
or two or three far-held measurements in the subsequent uniqueness study, 
and these are known to be the minimum number of measurements required. 

4.1. Uniqueness results. 

Definition 4.1. A bounded open subset of a two-dimensional plane is called 
a cell. A set D in is said to be a polyhedral scatterer if D is compact, 
M.^\D is connected and dD consists of hnitely many closures of cells. A 
polyhedral scatterer D is said to be a polyhedral obstacle if D is open. 

In what follows, we let and denote the classes of polyhedral scat¬ 
terers and obstacles, respectively. Clearly, ^s- We note for clarity 

that for a scatterer in it may contain a component which is solely a cell 
(it is referred to as a screen-component), whereas for an obstacle in all 
of its components must be solid. It is easily seen that any [D, III] or [D, IV] 
contained in is admissible which satishes the geometric conditions (12.71) 
and (|2.45p . If D S then the corresponding forward elastic scattering 
problem is well-understood, completely covered in our earlier discussion. 
However, if D £ we cannot hnd a convenient reference for the well- 
posedness of the corresponding forward elastic scattering problem. In the 
sequel, when a scatterer D £ t3^g is considered, we always assume the unique 
existence of an solution to the forward scattering problem, and 

focus on the study of the corresponding inverse problems, which is the main 
theme of the current section. Moreover, in order to ease the discussion, we 
present the uniqueness results for general polyhedral scatterers, whereas we 
only discuss the stability results for polyhedral obstacles. 

Theorem 4.2. Let D and D' he two polyhedral obstacles in such that 
U^ix] [D, K]) = U^{x; [7*^, [D', K']) Vx £ S^, (4.11) 

where t/*”’ is given in (|1.3p with either ap ^ 0 or ag ^ 0, and K, K' = HI or 
IV denote the respective kinds of the obstacles. Then one has 

[D,K] = [D’,K']. 


(4.12) 
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That is, both the obstacle and its physical type can be uniquely determined by 
a single total far-field measurement. Furthermore, if ap 0, in (14.111) 
can he replaced by U^; and if ag 0, in (14.lip can he replaced by Uf°, 
then the same conclusion holds true. 

Proof. First, we assume that ap 0. By the correspondences in Theorem l4.ll 
and Proposition 14.11 we have from (j4.11l) that 

Vf{x-[D, C(K)]) = V^{x- [D', C(K')]) Vx e (4.13) 

for K, K' = III or IV, with C(K) given by (14.3p . Then by the corresponding 
uniqueness result in US] for inverse acoustic scattering problems, one readily 
deduces (14.121) . Next, in the case ag 0, we have from (|4.11l) that 

Ef{x-ET, [I?,B(K)]) = Ef{x-ET, [I?',B(K0]) Vx G (4.14) 

where K,K' = III or IV, and the notation B(K) is defined in ()4.3p . Now, 
the uniqueness can be readily deduced from the corresponding results for 
inverse electromagnetic scattering problems in m- The other cases with 
Uf° replaced by U!^ or can be treated in a completely similar manner. 
The proof is complete. □ 

Theorem 4.3. Let [T*,K] and [D',K] be two polyhedral scatterers in l^g 
where K = III or IV represents the type of the scatterers. Then one has 

D = D\ (4.15) 

provided 

Uf°{x-,W^,[D,K]) = Uf°{x;U^^,[D',K]) Vx G (4.16) 

hold(s) for one of the following cases: 

(i) K = IV and ap 0, along with one single incident field C/*” = 
[/*”'(•; d, rf-*-, Op, Os, w) given by (II.3p .- 

(ii) K = IV and ag / 0, with two incident fields given by 17®” = 
W"'{-;d,dj-,ap,ag,uj), j = 1,2, such that d and dj-, j = 1,2, are 
linearly independent; 

(hi) K = III and ap 0, with the incident t/®” = W'^{-;dj,d-^,ap,ag,uj), 
j = 1,2,3, such that di, d 2 and ds are linearly independent; 

(iv) K = III and ag / 0, with t/®” = t/®®®(-; d, dj“, Op, Og, w), j = 1,2, 
such that di and d^ are linearly independent. 

Furthermore, if ap 0 (namely, in the cases (i) and (Hi)), in (I4.16p 
can be replaced by ; and if Og / 0 (in the cases (ii) and (iv)), Lff° in 
(14.lip can be replaced by , then the same conclusion holds true. 

Proof. Similar to the proof of Theorem 14.21 the primary strategy here is 
the application of Theorem 14.11 and Proposition 14.11 along with the use 
of the relevant uniqueness results for inverse acoustic and electromagnetic 
scattering problems. 

If K = IV and ap 0, similar to (I4.13p . we hrst obtain that 

v^{x-E^, [D,SS]) = v^{x-,v;^, [17',SS]) Vx G (4.17) 

with Up” given by (14.3p . By the corresponding uniqueness result in inverse 
acoustic scattering problems for sound-soft polyhedral scatters (cf. [XlfUji. 
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one readily has (14.151) . In the case that K = III and ap ^ 0, we obtain from 
(14.4p that 


i;“(x;i;;",[II,SH]) = i;“(x;n;^.,[Zl',SH]) Vx e S^, j = 1,2,3, (4.18) 


where 

:= -V • = -{iapkp)e^^^^-‘^\ j = 1,2, 3, (4.19) 

with three linearly independent directions di, d 2 and d^. Hence, (14.151) 
can be shown by the uniqueness results for inverse sound-hard problems 

(c.f. mi). 

If K = IV with Us 7 ^ 0, we have from Theorem 14.11 and Proposition 14.11 
that 


Ef{x-, Ef, [D, PMC]) = Ef^{x] Ef, [O', PMC]) Vx G S^, (4.20) 

with 

El^{x) = V A C/*"'(x; d, dj~, ctp, as, a;) = iaskg {d A dj~)e^^‘’^'^, j = 1, 2. 

(4.21) 

Define 

Es := -^V A Es(-; El\ [D, PMC]), K - "fv A Es(-; [D', PMC]). 

Kg f^s 

(4.22) 

Similar to (I2.32p - (|2.34l) . one has that {Eg, —Es{-; El"-, [D, PMC])) solves the 
Maxwell system (|2.25p with £'*"■ given in (I2.33p . and the PMC boundary con¬ 
dition in (|2.25p on dD replaced by the PEC one (|2.34p : and {E'^, — Es(-; E*"', [D', 
PMC])) solves the same problem except for the scatterer D replaced by D'. 
That is. 

Eg = Es(-;Ef, [D,PEC]), E' = Eg{-;Ef, [D',PEC]), (4.23) 

with two incident plane waves 

Ef = -iogkg j = 1,2. (4.24) 

Therefore, from (I4.20p . (I4.22p and (|4.23p . one observes that 

E^{x-,Ef, [D,PEC]) = Ef{x-,Ef, [D',PEC]) Vx G (4.25) 

for E*” given in (14.241) with two linearly independent vectors d^ and d^- 
Hence, (I4.15p can be obtained by the uniqueness results for inverse PEC 
problems in [ 12 | . 

In the last case for K = HI and as 7 ^ 0, similar to (I4.2np . one observes 
that 


Er(x;E-, [APEC]) = Er(x;E-,[Z)',PEC]) Vx G j = l,2, 

(4.26) 


with 


Ei^^{x) = iagkg {d A j = 1,2. (4.27) 

It is easy to see that d Adi and d A d 2 are linearly independent, since di and 
d 2 are linearly independent and both are perpendicular to d. Now, (14.151) 
can be deduced in the same way as that for the case (iii). 

Einally, for the other cases with replaced by or U^, one can show 
the corresponding uniqueness results by using completely similar arguments 
as above. The proof is complete. □ 
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We would like to point out that, the two uniqueness results in [5] in de¬ 
termining an elastic polyhedral scatterer are, respectively, case (i) and part 
of case (iv) in Theorem 14.21 A crucial ingredient in [5] to prove the cor¬ 
responding uniqueness results therein is the so called “reflection principle” 
for solutions to the Lame system dni), which was developed from the corre¬ 
sponding reflection principles for the Helmholtz equation and the Maxwell 
equations, respectively (cf. [UlISKTl]). By using the decoupling results in 
Theorem l4.ll the next proposition gives a simple proof of the reflection prin¬ 
ciple for the Lame system. It can be straightforwardly used to explain the 
reflection principles in [5] that were critical in the proofs of the relevant 
uniqueness results. Moreover, the reflection principle might find important 
applications in other context on elastic wave scattering. 


Proposition 4.2. Let the admissible scatterer D, and the functions U, Eg 
and Vp, he the same as in Theorem EH Let n be a two-dimensional plane, 
and let H' be the affine projection of H, namely, H' he the parallel plane of 
n passing through the origin. Set 


Lf (x) '.= RyiiU i^Rnx), x G RuD'', 


(4.28) 


where D'^ is given by := \ D, and Rjj represents the geometric reflec¬ 

tion with respect to the plane H. Then one has 


and 


with 


U — Up -\- Us 

Kp rvg 


Vp = —V • U and Eg :='V A U, 


(4.29) 

(4.30) 


Vp{x) := Vp{Ru{x)), Eg{x) := -Ru_'Eg{Rux), x G RnD'^. (4.31) 

Proof. It is verihed with straightforward computations that the reflected Vp 
defined in (|4.31l) is invariant under gradient, namely, 

Vvp{x) = Rw (Vup) (/?na;). (4.32) 

Similarly, Eg given by (|4.31l) is invariant under the curl operator, that is, 

V A Eg{x) = Ryi' (V A Eg) [Ryix). (4.33) 

Thus, by applying (14.2p one can derive for any x G RnD'^ that 


U{x) = Ryi'U{Rux) 

= 73 (V A Eg) (Rnx) (Vup) (i?na;) 

Kg Kp 


(4.34) 


Taking gradient and curl, respectively, on both sides of (|4.34p . one readily 
obtains (|4.3Up . 

The proof is complete. □ 
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4.2. Stability results. In this section, we present the stability estimates in 
determining polyhedral obstacles by a single far-held measurement, which 
quantihes the uniqueness results in Theorem 14.21 

First, we introduce some preliminary dehnitions and statements before 
presenting the main stability results. 

Given re G N, we say that a domain C M" is Lipschitz with constant r 
and L if the following assumption holds. For any x G 512, we hrst change 
rigidly the local coordinates such that x = 0. There exists a function ipx ; 

—>■ M, which is Lipschitz with Lipschitz constant bounded by L, such 
that (/9(0) = 0, and 

Brix) r\n = {y = {y', yn) G G M, y' G yn < 

and consequently, 

Br{x) ndQ. = {y = {y', y„) G Br{x);yn G M, y' G = ip{y')}. 

A scatterer C is said to be polyhedral with constants h and L, if 
dD is composed of a hnite union of the closures of cells, 

m 

dD=(jC'j, (4.35) 

i=i 

where rre G N is a certain hnite (unknown) number, and each Cj, 1 < j < rre, 
is a Lipschitz domain in with constants h and L; and two different cells 
may intersect only at boundary points. 

In this paper, the distance between two scatterers D and D' in is 
measured by the Hausdorff distance dn, 

dniD, D') := max < sup dist(x, D'), sup dist(x, D) > . (4.36) 

x£D' J 

Given positive constants r, L, R and h, let = ^g(r, L, R) denotes the 
set of obstacles D C Bji such that, 12 is a Lipschitz domain with constants 
r and L, and that D is polyhedral with constants h and L. In the rest of 
this section, we shall aways let the positive constants r, L and R, and the 
elements Up, Og, d, d^ and w mentioned in (II. 3p . and the wave numbers kp 
and kg be hxed; they are all referred to as the a priori data. We also keep 
the other positive number h being hxed. It is noted that h is not included 
into the a priori data. 

Theorem 4.4. Let [D, K] and [D', K] be two polyhedral obstacles in 3^^ = 
L, R), where K = III or IV represents the type of the scatterers. Sup¬ 
pose that 

[Z2,K]) - ^7^(•;G*^ [I2',K])||^2(s 2) < e, (4.37) 

for a single incident wave t/*"' given by (1131). Then there exists a positive 
number e{h;K), depending only on h, on the a priori data, and on the 
scatterer type K, such that if e < e{h;K), then one has 


dH{D,D')<C[iPise)r 


( 4 . 38 ) 
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where s is a certain number depending only on the a priori data which shall 
be specified below; the function if : {0,1/e) ^ (0,1) is defined by 

if{t) = exp(^-(log(-logt))^/^) ; 

C is a positive constant depending only on the a priori data and the scat- 
terer type K, but indpendent on h; and the positive constant a depends only 
on the a priori data and the scatterer type K, on or not on h; with more 
specifications as follows, 

(i) if K = IV and Up 0, then s = kp, the constant a in ()4.38l) can 

be chosen so that it depends only on the a priori data, not on h; 

moreover, in ()4.37p can be replaced by Ufi°; 

(ii) if K = III and Op / 0, then s = kp, the number a in (|4.38p does 

depend on h in general; and in (|4.37p can also be replaced by 

JJOO. 

Proof. Let be a radiating solution to the Lame system 0, and let 
and be the corresponding far-field pattern given by (|1.18l) and 
(I1.19p . First, by (ll.20p one observes that 

II^I°°IIl2 (§2 ) = ||L'^||j;^2 (§2 ) -I- ||L'“||l2 (S2). (4.39) 

It is clearly that Ufi-] C/*”, [D,K\) — Ufi-', [D' , K]) is a radiating solution 
to the Lame system (|l.ip in \{D D D'). Thus, (14.371) along with (j4.39l) 
implies that 

^(r,K) := ||[/r(-;C7*", [D,K]) - t/*", [77', K])||^2(s2) < e, (4.40) 

for T = p,s and t. 

Denote 

VpiK):=-V-U{-,W^,[D,K]) and v'p{K) := -V ■ ,[D',K]). 

(4.41) 

It is deduced by Theorem 14.11 that 

Vp{K) = vp(v^^; [D, C(K)]) and ^(K) = [D', C(K)]), (4.42) 

with u*** and C{K) given by (14.3p and (14.51) . respectively. Therefore, one 
observes from Proposition 14.11 that 

||u“(K) - u;°°(K)||i2(s2) < kpap,K) < kpE, (4.43) 

where u))°(K) Up°°(K) are defined similarly to v’^ in Proposition 14.11 Now, 
the stability results in case (i) can be obtained by the corresponding stabil¬ 
ity estimate of inverse acoustic scattering for polyhedral sound-soft obstacles 
in [20]; and in case (ii) can be concluded by the stability estimate for deter¬ 
mining sound-hard obstacles in HU. 

The proof is complete. 

□ 


5. Concluding remarks 

In this paper, we consider time-harmonic elastic wave scattering governed 
by the Lame system associated with a third or fourth kind scatterer. It is 
well know that the elastic wave field can be decomposed into the pressure 
(longitudinal) and shear (transversal) parts. Generally, the pressure and 
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shear waves coexist, but propagating at different speeds. We derive two 
geometric conditions, related to the mean and Gaussian curvatures of the 
boundary surface of the underlying scatterer, which can ensure the com¬ 
pletely decoupling of the pressure and shear waves. Dehnitely, this phenom¬ 
enon is of signihcant physical interest. On the other hand, if the decoupling 
occurs, then mathematical reduction occurs where the elastic wave equation 
can be decoupled into two simpler Helmholtz and Maxwell systems. Us¬ 
ing the decoupling results, we obtain comprehensive uniqueness results for 
inverse elastic scattering problem in determining polyhedral scatterers by 
the minimal number of far-held measurements. Furthermore, we quantify 
the uniqueness result in determining polyhedral obstacles to derive opti¬ 
mal stability estimate of logarithmic type by a single far-held measurement. 
This is the hrst stability result in the literature for inverse elastic obstacle 
scattering problem by a single far-held measurement. We believe that the 
decoupling results can hnd more important applications in other context. 
For example, in [8], the invisibility cloaking of elastic waves was considered. 
Since the pressure and shear waves propagate at different speeds, one may 
hrst decouple them and then cloak them separately. In doing so, more fa¬ 
vorable cloaking devices can be achieved. We shall explore such an idea in 
a forthcoming article. 
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